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Abstract 

The neutron capture cross section of the unstable nucleus 186 Re is studied by investigating the 
inverse photodisintegration reaction 187 Re( 7 ,n). The special interest of the s-process branching 
point 186 Re is related to the question of possible s-process contributions to the abundance of the 
r-process chronometer nucleus 18 'Re. We use the photoactivation technique to measure photodis¬ 
integration rates. Our experimental results are in good agreement with two different statistical 
model calculations. Although the cross sections predicted by both models for the inverse reaction 
186 Re(n, 7 ) is too low to remove the overproduction of 186 0s; the two predicted neutron-capture 
cross sections differ by a factor of 2.4; this calls for future theoretical study. 
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I. INTRODUCTION 


Q- 


Almost all elements above mass A « 60 can be produced in neutron capture reactions 
Two different neutron induced processes are necessary to explain the abundance distribution 
of heavy elements. The first one is the slow neutron capture process (s-process). The neutron 
densities are of the order of n n & 10 8 cm -3 and the time scale r n between two subsequent 
neutron capture reactions is typically of the order of years. The s-process path propagates 
along the valley of stability. Whenever an unstable nucleus with a mean lifetime r < r n is 
reached, this nucleus /9-decays. If r ~ r n , a branching occurs and the s-process path splits. 
Thus, nuclei with r ~ r n are called branching points of the s-process. The second process 
is the rapid neutron capture process (r-process). High neutron densities ( n n 10 2O cm -3 ) 
lead to the production of very neutron rich nuclei up to 20 mass units away from stable 
nuclei. During freeze out, these nuclei /3-decay back to the valley of stability. 

There are at least two scenarios known where the s-process takes place. It occurs during 
helium burning in red giant stars and during helium shell flashes in low mass asymptotic 


giant branch (AGB) stars Q . The former scenario is mainly responsible for the production 
of elements between iron and yttrium. The latter, for the production of elements between 
zirconium and bismuth. For a detailed discussion see e.g. J]. In the following we will focus on 
the mass region A & 185 and, hence, restrict our discussion to the so-called main component 
of the s-process. 



V y V V \ r-process 


FIG. 1: The s-process path in the W-Re-Os mass region. Unstable nuclei are marked by dashed 
boxes (except 18 'Re). The indicated values are laboratory half-lives. However, the half-life of 18 'Re 
decreases by 10 orders of magnitude at typical s-process temperatures of T = 3 x 10 8 K and 18 'Os 
becomes unstable 


Due to its very long half-life (ti/ 2 = 5 ■ 10 10 a) the nucleus 18 'Re can be used as a r-process 
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chronometer |8|, |9| . The ratio iV( 18 'Re)/iV c ( 18 'Os) is related to the starting point of the r- 
process in our galaxy and, hence, to its age. N denotes the total and N c the cosmoradiogenic 
part of the abundance stemming from the decay of 18 'Re. To extract the cosmoradiogenic 
part of the 18 'Os abundance one has to subtract the s-process abundance N s from the total 
abundance N. In Fig. El the s-process flow through the W-Re-Os isotopes is shown. The 
s-process abundance of 18 'Os can be derived from the abundance of the neighboring s-only 
nucleus 186 0s via the local approximation jlf)| : 


iV s ( 187 Os)/iV s ( 18fi Os) « Fb n ( 18e Os)/b n ( 187 Os), 


/186 


- /' 186 / 


I- /187 f 


( 1 ) 


where <r n are the Maxwcllian-averaged radiative neutron capture cross sections (MACS) 
from the ground state, and F accounts for the correction of the cross section due to neutron 
capture on thermally excited states in 18 'Os, in particular on the first excited state at 
9.75 keV. TMs correction factor was flrSt calculated Q (see Q] for discussion). The 
neutron capture cross sections of 186 Os and 18 'Os were measured by Browne & Berman 
Browne, Larnaze & Schroder [kJ. by Browne & Berman Q] and Winters & Macklin [l5j]. 
resulting in an uncertainty of about 20% for the ratio R = CT n ( 186 0s]/<7 n ( 187 0s). Recently, 
these cross sections were measured by the n_TOF collaboration [la, Il7 |. 

The use of the Re/Os clock is not free of problems. First of all, the half-life of 187 Re 


H, 


strongly depends on temperature and 18 '0s becomes unstable under stellar conditions 
0 7|. Thus, it is necessary to use chemical evolution models of the galaxy f| to include 
irradiation effects on the abundance ratio R. The two branchings at 185 W and 186 Re also 
affect the s-process abundances in this region. Finally, the N s a n correlation Eq. m for the 
two s-only isotopes 186 ’ 18 'Os is not fulfilled. This can be caused by two facts. Either the 
branchings are not correctly modeled or the capture cross sections are strongly affected by 
stellar conditions. The branching at 185 W has already been studied and an overproduction 
of 186 0s was reported due to the new experimental value [3]. Thus the radiative neutron 
capture cross section of 186 Re is the only relevant cross section in this mass region which is 
not known experimentally yet. 

In this paper we study the radiative neutron capture cross section of the branching point 
nucleus 186 Re using an indirect method. The unstable nucleus 186 Re decays via (3 ~-decay to 
186 0s or via electron capture to 186 W with a half-life of t\/ 2 = 3.7d. Due to the fact that 
neutron capture experiments with such short-lived targets are nearly impossible, we choose 
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the inverse reaction 18 'Re( 7 ,n) 186 Re for onr investigation. After neutron emission the 186 Re 
nucleus is in the ground state or some excited state and the measured cross section is a sum 
over several channels. With maximum excitation energies just above the neutron threshold, 
only the lowest states can be reached, e.g., the first excited state in 186 Re at 59 keV. On 
the other hand, these low lying states in the 186 Re nucleus are also thermally populated 
- however not in the same proportions - under s-process conditions and contribute to the 
neutron capture cross section. Thus, the 18 'Re( 7 ,n) 186 Re cross section and the cross section 
of the inverse reaction 186 Re(n, 7 ) 18 'Re are related via the principle of detailed balance. 

In section QT] we describe our experimental setup. Section |HT| explains the analysis of our 


data and the results for the 187 Re( 7 ,n) cross section are presen 


to calculations using the NON-SMOKER Jl9j, |20j and MOST 


;ec 


21 


. The results are compared 


Eg 


22] codes. Both computer 


codes are based on the statistical Hauser-Feshbach model but use different input parameters. 
In section |Tvl the implications on the Re/Os clock are briefly discussed. 


II. EXPERIMENTAL SETUP 



FIG. 2: The real photon setup at the superconducting electron accelerator S-DALINAC at Darm¬ 
stadt University of Technology. The electron beam with a maximum energy of 10 MeV and maxi¬ 
mum intensity of 40 /xA hits a thick copper target and produces bremsstrahlung. The beam position 
is monitored with an ionization chamber. The flux is monitored online with two actively shielded 
100% HPGe detectors using the n B( 7 , 7 ’) reaction. 
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The 18 'Re( 7 ,n) 186 Re experiment was performed using the photoactivation technique at 


the real photon setup 


23] at the superconducting electron accelerator S-DALINAC |2J] at 


Darmstadt University of Technology (see Fig. The monoenergetic electron beam hits 
a thick copper radiator (d ~ 1.4 cm), where it is completely stopped and converted into a 
continuous spectrum of bremsstrahlung photons with an endpoint energy E max . The photons 
are collimated and irradiate the target of interest at a distance of d ~ 150 cm behind the 
radiator. This leads to a well-defined photon beam with a spectral composition that is 
analyzed in detail, see e.g. j3]. The y-intensity as well as the electron current are monitored 
online in order to control the beam position on the copper radiator. The targets consist 
of thin metallic rhenium discs (m ~ 340 mg, 0 = 2 cm) of natural isotopic composition 
and of two layers of boron (m ~ 650 mg each, 0 = 2 cm) embedding the rhenium disc 


with a sandwich-like structure. The incoming p 
nuclear resonance fluorescence reaction 11 B( 7 , 7 ’) 


roton intensity is normalized using the 
3] (see Fig. EH). The scattered photons 


are registered online with two HPGe detectors positioned at 90° and 130° with respect to 


the beam direction. The detector efficiency e was measured up to 3.6 


calibration sources. We have used the Monte Carlo code GEANT 
efficiencies up to an energy of 10 MeV. 


VleV using standard 


27] to extrapolate the 



FIG. 3: Adjustment of the simulated spectral photon distribution to the intensities IV 7 obtained 
from the reaction 11 B( 7 , 7 ’) at E max = 9.9 MeV. The flux is integrated over 23 hours, a typical 
activation time. 


The spectral photon distribution was simulated with the same code and is adjusted 
to the photon intensities obtained from the n B( 7 , 7 ’) reaction. The result for the case 
E max = 9.9 MeV is shown in Fig. 01 A list of all measured energies is shown in Tab. ED The 
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E 7 [keV] 


FIG. 4: Typical 7 -spectrum measured after photoactivation of nat Re. The line at 122 keV and at 
137keV stem from the /3-decay of 186 Re into 186 W and 186 0s, respectively. The line at lllkeV 
stems from the decay of 184 Re into 184 W. The target was activated for 23 hours at an energy of 
£max=9.9MeV and the 7 -spectrum was accumulated over a period of 23 hours starting one hour 
after activation. 


statistic uncertainty for the photon flux calibration is very small because of the well known 
photo-response of n B Q], The systematic uncertainty includes a 5% uncertainty from the 
detector efficiency extrapolation to higher energies and the uncertainty dne to the shape of 
the photon spectrum close to the endpoint energy. This uncertainty is smaller for higher 
endpoint energies because the 7 -transitions at 8.92 MeV in the reaction 11 B( 7 , 7 ’) can be 
used for calibration if F/ max > 9 MeV. The total uncertainty for the photon flux calibration 
can be estimated to be about 12 % to 21 % depending on the endpoint energy E max (see 

Tab.HH). 


After activation, the 7 -rays from the decay of the unstable nucleus produced are measured 
offline using a well shielded HPGe detector with an energy resolution better than 0.15 % and 
an efficiency of 30 % relative to a 3" x 3" Nal detector. The detector efficiency is determined 
between 60keV and 1.3 MeV with standard calibration sources. The target geometry and 


self absorption effects are simulated with a Monte Carlo code |27[. A typical spectrum after 
activation is shown in Fig. 01 

The peak area Y is directly proportional to the number of ( 7 ,n)-reactions TZ in the target 
during the activation. Knowing the detector efficiency e and the absolute 7 -branching of 
the transition, the factor of proportionality can be directly calculated for each run using the 


6 







TABLE I: The peak area Y of the two 7 -transitions following the /3-decay of 186 Re and the 
corresponding number of ( 7 ,n)-reactions 77. The results for the transition at 122 keV at E max = 
7.65 MeV are omitted due to low statistics. 


Bmax [MeV] 

Y(122 keV) [10 3 ] 

Y(137keV) [10 3 ] 

77(122 keV) [10 6 ] 

77(137 keV) [10 6 ] 

9.9 

35.12 ±0.22 

605.7 ± 0.8 

287 ± 14 

286 ± 14 

9.45 

17.48 ±0.15 

301.3 ±0.6 

148 ±8 

147 ±7 

9.0 

9.91 ±0.20 

167.6 ±0.4 

85.0 ±4.6 

83.0 ±4.2 

8.55 

3.36 ± 0.07 

58.81 ±0.25 

33.7 ±1.8 

34.1 ± 1.7 

8.325 

1.79 ±0.05 

29.66 ±0.18 

16.3 ± 1.0 

15.6 ±0.8 

8.1 

0.87 ±0.04 

14.98 ±0.18 

6.86 ± 0.48 

6.85 ± 0.35 

7.875 

0.50 ± 0.04 

7.62 ± 0.09 

3.24 ±0.31 

2.85 ±0.15 

7.65 


6.82 ±0.08 


0.39 ± 0.02 


law of exponential decay. The number of ( 7 ,n)-reactions 77 is proportional to the energy- 
integrated cross section I a , 

/“^max 

77 = n t I a = n t / N^E, E max )a(E) dE (2) 

where n t is the number of 187 Re atoms, S n is the neutron threshold, E max is the energy 
of the electron beam and, hence, the endpoint energy of the spectral density distribution 
iV 7 (£’, _E max ) of the bremsstrahlung photons, and cr(E ) is the 18 'Re( 7 ,n) 186 Re cross section. 
We used the strong transitions at 137 and 122 keV, respectively, for our analysis. 

The uncertainty of I a is about 5% to 6 %. The statistical uncertainty is small due to 
the high sensitivity of the photoactivation technique (see Fig. ED. The uncertainty of I a is 
dominated by the 5 % uncertainty stemming from the simulation of the detector efficiency 
including self-absorption effects in the target. The peak areas Y and the number of ( 7 , 11 )- 
reactions 7 Z are summarized in Tab. [I] The different components of the uncertainties as well 
as the total uncertainty are summarized in Tab.lIIlfor each endpoint energy 77 max . 


III. ANALYSIS OF THE DATA 

The peak area in the activation spectrum is only proportional to I a , thus, we cannot 
determine the cross section directly. Assuming a certain shape of the cross section, one can 
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TABLE II: Summary of experimental uncertainties for the different endpoint energies E max . The 
details are discussed in the text. 


E/max [keV] 

7 -intensity 

stat. [%] 

sys. [%] 

shape[%] 

number of ( 7 ,n)-reactions 

stat. [%] sys. [%] 

Total [%] 

9900 

1.0 

9.8 

8 

0.1 

5 

13.6 

9450 

1.4 

8.6 

8 

0.2 

5 

12.8 

9000 

1.2 

10.7 

11 

0.2 

5 

16.2 

8550 

1.3 

12.9 

11 

0.4 

5 

17.7 

8325 

1.5 

10.7 

11 

0.6 

5 

16.2 

8100 

2.0 

14.1 

14 

0.8 

5 

20.6 

7875 

1.8 

14.9 

14 

1.2 

5 

21.2 

7650 

2.0 

13.0 

14 

1.2 

5 

19.9 


verify and normalize this assumption by measuring at several endpoint energies because the 
centroid of the integrand in Eq. © changes for each endpoint energy. 

In case of pure s-wave neutron emission, the ( 7 , 11 ) cross section can be parametrized 


as 


29] 


cr(E) = cr 0 


E-Sr, 


S n 


(3) 


for the energy region close above the neutron threshold energy S' n ( 187 Re) = 7.363 MeV. The 
combination of Eq. 0 with Eq. 0 yields the normalization factor 07 . The normalization 
factors shown in Fig.lHl as a function of endpoint energy should be constant. We have 
measured the integrated cross section using bremsstrahlung with endpoint energies of 7.65, 
7.875, 8.1, 8.325, and 8.55 MeV, respectively. For S n < E < 8.55 MeV we obtain: 


a(E) 


(80.4 ± 9.6) mb 


E S n 

Sn 


(4) 


Between 9.5 and 19 MeV, the ( 7 , 11 ) cross section can be parametrized by a superposition 
of two Lorentzians Q. Between 8,5 and 9 ,MeV, we interpolate tire cross section with a 
third order polynomial. The requirement that the cross section as well as its first derivative 
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are continuous, determines the four parameters of the polynomial: 


°{E) 

mb 


= -17 


- 4.3 x 10 


E \ 3 
MeV ) 

3 / E ^ 


+ 4.7 x 10 2 


f—y 2 

VMeV J 


+1.3 x nr 


( 5 ) 


MeV J 

The parametrization of the cross section for the energy range between S n and 10.5 MeV is 
shown in Fig.|Hl The measurements performed at E max = 9.0, 9.45 and 9.9 MeV verify this 
parametrization up to the low energy tail of the giant dipole resonance. Our data fit nicely 


with an older experiment jsflj], in which the 18 'Re( 7 ,n) cross section was measured at higher 
energies around the giant dipole resonance. 

Alternatively one can start with a cross section er t h calculated with the NON-SMOKER 
or MOST code. A normalization factor / is introduced for absolute calibration: 


I* = f' 


a th (E)NJE,E n 


dE 


( 6 ) 


The normalization factor / needs to be energy independent presuming that the predicted 
shape of the cross section cr t h is correct. Additionally, if the absolute value of the calculated 
cross section is correct, / should be close to 1.0. The normalization factors / are determined 
to be 


/nons = 1.15 + 0.31 

( 7 ) 

/most = 1.07 + 0.28 

(8) 


The normalization factors as a function of the endpoint energy E max are shown in Fig.0 
and the normalized cross sections are shown in Fig.|Hl The agreement between the two 
calculations and our experimental cross section is remarkably good. The small deviations 
in the threshold region result either from experimental uncertainties regarding the precise 
shape of the bremsstrahlung close to the endpoint energy E max or from numerical problems 
occurring in the theoretical calculation. 


IV. IMPLICATIONS FOR THE RE/OS CLOCK 


The current uncertainty of the time duration of nucleosynthesis obtained from the Re/Os 


clock is about 2.3 Gyr 


31]. The main source of uncertainty with regard to nuclear physics 


aspects of the clock are presently the neutron capture cross sections. 
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FIG. 5: The normalization factors co of Eq. ©, alternatively / of Eq. © are shown as a function 
of the endpoint energy F? max of the photon spectrum for our experimental cross section and for the 
two theoretical predicted cross sections. Both model predictions seem to underestimate the cross 
section slightly in the vicinity of the reaction threshold. 

The Maxwellian-averaged neutron capture cross sections in Tab. lllll are derived using the 
MOST and NON-SMOKER code. The ( 7 ,n) cross sections are in good agreement with the 
present experimental results (see Fig. ED- However, the neutron capture cross sections of the 
inverse reaction differ by a factor of 2.4. 

Taking the value —ended by Bao el Q ike branehing ratio between .-deeay 
and neutron capture is R — \ n /\p = 5.4 x ICE 4 . Considering the values from the NON- 
SMOKER and MOST calculations, one obtains Rnons = 5.2 x 10~ 4 and Rmost = 2.1 x 


i 1 i ' i ' i 

& experimental cross section 


0 NON SMOKER- 


<> ° 0 


0 <> 


<> 0 0 


0 MOST- 


J_i_L 
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FIG. 6: Comparison of the experimental 187 Re( 7 ,n) cross section to the calculated cross sections 
from the MOST and NON-SMOKER code. Additionally, the data points extracted from Goryachev 


et al. [3d] are plotted. The shaded area represents the uncertainty of our experimental results. 


TABLE III: The MACS calculated with an updated version of NON-SMOKER [20] and MOST 
for a typical s-process temperature of kT = 30 keV. The stellar enhancement factor as well as the 


adopted values of the MACS were taken from Bao et al. 


321 



ctgs [mb] 

^thermal [mb] 

NON-SMOKER [2Q] 

1485 

1546 

MOST [21] 

616 

623 

Ref. [32] 

1550 ± 250 

1615 ± 260 


10 , respectively. Both models predict a (n, 7 )-cross section which is smaller than the 


value recommended 


in 3- 


This would lead to a further enhanced production of 186 0s. 
Nevertheless due to the small branching the abundance distribution of the elements in this 
mass region is not changed significantly. 


V. SUMMARY AND CONCLUSION 

We have measured the 18 'Re( 7 ,n) cross section just above the reaction threshold. The 
cross sections calculated within the statistical model using the computer codes NON- 
SMOKER and MOST are in good agreement with our experimental data. Even though 
the predictions for the photodisintegration cross sections are nearly identical, the MACS for 
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the neutron capture reactions differ by a factor of 2.4. At this point, further theoretical 
investigations are necessary. 

The MACS calculated with the updated NON-SMOKER code is close to the value given 
in ( 3 ^. Thus, the overproduction of 186 0s reported in remains. The MACS calculated 
with the MOST code is significantly lower and, therefore, a further increase of the 186 0s 
production is predicted. For a more quantitative statement, complex network calculations 
are mandatory. 

The disagreement between the two models may have its origin in the lack of precise 
nuclear data that enter into both models. This shows again the need for high precision mass 
measurements, El-strength distribution studies, and the study of nucleon-nucleus optical 
potentials. 

The fact that both model predictions are in good agreement with our data but could not 
reduce the overproduction of 186 0s supports the idea that the adopted value of the 186 0 s(n, 7 ) 


cross section 


ity at CERN 16 


32] is too small. Very recent results from an experiment at the n-TOF facil- 


and at the Forschungszentrum Karlsruhe [IjJ point in this direction as well. 
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